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ABSTRACT 



The intent of this paper is to describe a method of coloring a 
map and to present an algorithm for the solution of this problem . 

A computer program was developed to provide solutions to the problem 
of coloring a map which consists of a finite number of areas. This 
algorithm may also be applied to problems other than map-coloring. 
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I . INTRODUCTION 



When coloring a geographical map, it is customary to use different 
colors for adjacent regions - where adjacent means a ccmmon finite 
boundary, it has been conjectured, but not proved, that "For any sub- 
division of the plane into non-overlapping regions, it is always pos- 
sible to mark the regions with one of the numbers 0, 1, 2, 3, in such 
a way that no two adjacent regions receive the same number."^ 

The statement that a map can be colored with four colors is 
accredited to Moebius, who first proposed it in 1840. This conjecture 
has never been proven; however, a satisfactory counter-example has not 
been demonstrated either. A proof for five-coloring a map may be found 
in [Courant and Rcbbins 1941] and [Oystein 1967] . 

This thesis will use the four-coloring conjecture without proof 
and present an algorithm that utilizes an iterative routine to determine 
what color each region should be. A similar iterative technique for 
general integer programming problems is demonstrated in [Greenberg April 
1969] and [Greenberg May 1969], which present the proofs. We are 
interested in showing whether the iterative technique can be used for 
large-scale problems such as the map-coloring problem. 

The problem formulation suggested by Gcmory [Dantzig 1960] is used. 
Let the map regions be denoted by i = 1,2, ...,N and let x^ be an integer 
valued variable such that 

0 4 x i < 3 (1) 

^Courant, Richard and Rcbbins, Herbert; What is Mathematics? , p. 247, 
Oxford University Press, 1941. 
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where the four values x^ = 0, 1, 2, or 3 correspond to four different 
colors. Since any two adjacent areas must have different colors, the 
constraint may be written in an either/or form 



either x. - x. >, 1 
i D 



or x. - x. > 1 
D i 



(2) 



if i and j have a common border. Equation (2) may be rewritten as: 



x. - x. > 1 - 46. . (6. . = 0,1) 

i D 13 ID 



x . - x. >, -3 + 46 . . 
3 1 13 



(3) 



Instead of using (2) and (3) we can simplify the problem by using 



x. - x. ^ 0 
1 3 



(4) 



where i and j have a cannon border. 

This formulation can be written in matrix form 

Ax ^ 0 (5) 

where the elanents of a j_j' are either 0, -1, or 1. We can now 

enumerate all values of the left side of (5) . The column with the 

_ 2 

least number of zeroes is new selected from A^^^, , and the value 

of the corresponding x^ is set equal to one. The column a^ is then 
added column by column to the matrix to form a new matrix. The 

process continues in the same way by searching the new matrix for the 
column with the least number of zeros. In each iteration the correspond- 
ing Xj is increased by one. The column values represent the value of 
the left side of (5) . A solution occurs when a column is achieved that 
has no zeros. 



2 

The speed with which the algorithm reaches a solution is increased 
by modifying the method of choosing a . . This modification is explained 
in Section II. A. 4. 
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The algorithm presented is based upon the four-color conjecture, 
but is not restricted to the use of only four colors, if more colors 
are desired. 



II. THE COLORING PROCESS 



A. THE ALGORITHM 



1. Let x. = 0 for i = 1,...,N. 

i 

2. For every 2 regions of the map with a ccmmon border it is 
required that x^ - x^ ^ 0, i,j, = 1,...,N; i ^ j. 

3. Form the M x N matrix defined in step 2, where M is 

the number of equations required to describe the border constraints 
and N is the number of areas to be colored. The elements of ^MxN^ij' 
are either -1, 0, or 1. Note : 

a. For each row, a^, i = 1,. • . ,M, there are only two non zero 
elements, -1, and 1, corresponding to the constraints in step 2. 

b. For column a ^ , j = 1,...,N, has at least one non zero 
element and fewer than M non zero elements. 

4. Let ctj, j = 1,...,N, be the number of zero elements of each 

a^. Let 3 j , j — 1,...,N, be the number of zero elements of each a^ 
which appear in those positions where a^ = 0. Note : 3j < a ^ , 

j = 1, . . . ,N. Find the column of A^^ with the smallest 3 j . Let j** 

be the smallest j such that 3-** = min {3 }. Let j* = j**. 

3 j : 

5. Let x' ... = x. + + 1. Check all the constraints in which x' 

3* 3 * 3 * 

appears to insure that the constraint equations (4) are satisfied. If 

the constraint equations (4) are not satisfied, continue increasing 

x . . in increments of one until either (4) is satisfied or x' = 3. 

3 3 

If x' = 3 and (4) is not satisfied, set x' + = 0 and let j* be the 
: 3 



li 



smallest j > j** such that g.* = min {g.} and repeat step 5. If no 

3 j 3 

such j* exists, then find the smallest cu and proceed as with the 
3j's. If no ou exists such that x^* and x^ satisfy (4) for all i ^ j*, 
go to step 9. 

6. Determine the elements of column a^* by substituting the presently 
defined values of x^, i = 1,...,N, into the constraint equations (4). 



Add the column a^*, j = 1,...,N, to each column of A^^. Let A' 
^MxN + ^ a i*' * * * ' a j*J * 



MxN 



7. Let a., j = 1,...,N, be the number of zero elements of each a' .. 

3 J J 

Let 3^, j = 1,...,N, be the number of zero elements of each aV which 

appear in those positions where a^ = 0. Let j** be the smallest j such 

that 3. = min {g.} . Let j* = j**. 

3 3 3 

8. Continue steps 5 through 8 until cu =0, in which case a solution 

has been reached. Note : If x^ = 3 at any point during steps 5 through 

8, disregard the corresponding column a^ in any future iteration, since 
the corresponding x^ is at its maximum possible value. 

9. Let the column a. + be the one used to generate the latest A',, „ 

j* 3 MxN 

matrix. For each zero element of a^* there are 2 regions, g and h, with 
a ccutroon border where x^ - x^ = 0, i.e. , a violation of (4) . Let F be 
the set of all such regions g and h. Denote by the subscript f the 
element of F. 

10. For each f , increment the x^'s such that x^ = x^ + 1. Check all 
of the constraints in wliich x^ appears to see if the constraint equations 
(4) are satisfied. If for any f, an elonent of F, the constraint equations 
(4) are satisfied set f = j* and go to step 6. If (4) is not satisfied 
for any i, not an element of F, where regions i and f have a ccrrmon 
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Initial 




Possible Path to a solution to the map-coloring problem 



Figure 1 
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border, note the iteration that defined the present value of x^. 
Repeat step 10 until x^ = 3. Reset x^ = 0. Repeat step 10 for all 



11. Order the iterations defined in step 10 in decreasing 
iteration number. Let these iterations be called possible iterations. 

12. Reconstruct the most recent possible iteration and examine 
3^ = or ou = ou*, j = j* + 1,...,N, depending on which was used 
to determine x^* in the iteration being examined. Choose a new j* 

in the same manner as the j* was chosen in step 5. When a j* is found 
such that Xj* and x^ satisfy (4) for all i ^ j*, go to step 6. If no 
such j* exists satisfying (4) , then the iteration being examined is not 
possible. Repeat step 12. 

This canpletes the algorithm. If no possible solution exists, the 
algorithm will continue "backing up" the generated matrix chain until 
the original matrix is reached. Figure 1 shows a possible path to the 
solution of the map-coloring problem. 



B. EXAMPLE 




Fig. 2. Example Problem 



J The speed with which the algorithm reaches a solution is increased 
by starting the backtracking procedure when no variable associated with 
the minimum 3 . or a . can be set at any integer value between 1 and 3 . 
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where 



X 1 - x 2 * 0 
X 1 - x 3 * 0 
X 1 ' x 4 ^ 0 



X 2 “ x 3 ^ 0 
X 2 “ x 4 ^ 0 
x 3 " x 4 * 0 



set x. = 0 for all i 

l 

we write 



D 

5 . 

D 



X 1 


X 2 


X 3 


X 4 


1 


-1 


0 


0 


1 


0 


-1 


0 


1 


0 


0 


-1 


0 


1 


-1 


0 


0 


1 


0 


-1 


0 


0 


1 


-1 


3 


3 


3 


3 


3 


3 


3 


3 



column used 



then x^ = 1 

the column to add (a^*) is defined by: 



X-i 


- X„ = 


1 


x„ 


- x_ 


= 0 


1 


2 




2 


3 


X 1 


_ X 3 = 


1 


X 2 


" X 4 


= 0 


i — 1 

X 


1 

X 

(1 


1 


X 3 


" X 4 


o 

II 



where a^ = [111000] 



then adding [a^,...,! ] + we write 
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X 1 x 2 x 3 x 4 

2 0 11 

2 10 1 

2 110 

0 1-10 

0 10-1 

0 0 1-1 

a j 3 2 2 2 

Bj 3 1 1 1 



column used 



then = 2 

the column to add (a^*) is defined by: 



where a . ■ 

3 * 

then adding 



column used 



X 1 


X 2 


-1 


X, 


x n 


- X-, = 


1 


X, 


1 


3 






X 1 


i 

X 

li 


1 


X 



[- 111220 ] 

[a. 

*/•••/ ^ j * 1 + *MxN We write 



X 1 


X 1 


X 3 


X 4 


0 


-2 


-1 


-1 


2 


1 


0 


1 


2 


1 


1 


0 


2 


3 


1 


2 


2 


3 


2 


1 


0 


0 


1 


-1 


2 


1 


1 


1 


1 


1 


0 


0 


* 



X- 



X, 



X, 



= 2 
= 2 
= 0 
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then = 3 

the column to add (a.*) is defined by: 

3 



X-I 


- x n = 


-1 


x_ 


- x_, = 


-1 


1 


2 




2 


3 




x. 


- x_ = 


-2 


x_ 


- x, = 


2 


1 


3 




2 


4 




X 1 


ll 

X 

1 


1 


X 3 


I 

X 

II 


3 



where a.* T = [-1 -21-12 3] 

3 

then adding [a^*, . . . /a^] + we write 



X 1 x 2 X 3 X 4 

0 -2 -1 

-1 -2 -3 

2 11 
-1 0 -2 

2 3 2 

3 3 4 



a. 1 10 

3 

* 



there is no need to calculate column a' ^ because there are no zeros 
in column a' ^ and a solution has been reached. The solution is: 



Area number 



Color number 



1 

2 

3 

4 



1 

2 

3 

0 
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III. CONCLUSION 



The algorithm described will solve the map- coloring problem. 
Although only four colors were used in the example in Section II-B, 
the algorithm can be adapted to solve a problem for any feasible 
number of colors. Several representative map-coloring problems were 
solved using the enclosed computer program. Appendix A shows a 
solution to four- coloring the continental United States. In the 
representative problems solved, both four and five colors were used. 

If the solution is not feasible, e.g., the use of four colors 
when five colors are required, as in the split- state problem, the 
algorithm may consume several hours of computer time before showing 
that no solution exists. The computer time is consumed in back- 
tracking the solution path and in investigating every possible 
solution path for both the a ' s and the 3 ' s . 
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APPENDIX A 



SAMPLE PROBLEM SOLUTION 



The solution for four -coloring the continental United States of 
America , including boundaries and major water bodies, is as follows: 



STATE COLOR 

Maine 0 
New Hampshire 3 
Vermont 0 
Massachusetts 2 
Rhode Island 3 
Connecticut 0 
New York 3 
New Jersey 0 
Pennsylvani a 2 
Delaware 3 
Maryland 0 
Virginia 2 
North Carolina 3 
South Carolina 0 
Georgia 2 
Florida 0 
Alabama 3 
Tennessee 1 
West Virginia 1 
Ohio 0 
Indiana 2 
Kentucky 3 
Illinois 0 
Wisconsin 2 
Michigan 3 
Minnesota 0 
Iowa 3 
Missouri 2 
Arkansas 3 
Mississippi 2 
Louisiana 0 
Texas 2 
Oklahoma 0 
Kansas 3 
Nebraska 0 
South Dakota 1 
North Dakota 3 
Montana 0 
Wyoming 2 
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Colorado 1 

New Mexico 3 

Arizona 1 

Utah 3 

Idaho 1 

Washington 0 

Oregon 2 

Nevada 0 

California 3 

BORDER AREAS 

Canada 2 

Mexico 0 

Water - including: 1 

Atlantic Ocean 
Pacific Ocean 
Great Lakes 
Great Salt Lake 



The above problem was solved in 97.2 seconds on the IBM 360/67 computer 
and used 81 iterations. 
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APPENDIX B 



PROGRAM FLOW CHART 
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ITERATIVE ROUTINE TO FIND A SOLUTION TO THE MAP COLORING PROBLEM 
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